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ABSTRACT 

The Eulerian and Lagrangian second-order perturbation theories are solved explicitly in 
closed forms m ^ 1 and A 7^ Friedmann-Lemaitre models. I explicitly write the second- 
order theories in terms of closed one-dimensional integrals. In cosmologically interested cases 
(A = or i7 + A = 1), they reduce to elementary functions or hypergeometric functions. For 
arbitrary and A, I present accurate fitting formula which are sufficient in practice for the 
observational cosmology. It is reconfirmed for generic VL and A of interest that second-order 
effect only weakly depends on these parameters. 
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The gravitational instability presumably plays an essential role in the formation of the 
large-scale structure of the universe. The dynamics of such interesting phenomena involves 
the nonlinearity which is difficult to deal with. Any exact solution for nonlinear evolution 
is not known in general situation. We have been mainly resorted to A^-body methods for 
fully nonlinear problems in this field. Although such methods can shed light on strongly 
nonlinear regime, the resolution is fairly limited and they can barely survey only small fraction 
of parameter space of possible models. On large-scales where density fiuctuations are small, 
perturbation theories are quite powerful. They can analytically give results for large parameter 
space of possible models. 

Two formulations for higher-order perturbation theories are focused in the literature. One 
is Eulerian formulation |l|,||,||,H,|],|^,|^,||,||,|lO|,|ll|,|l^,|l3l and the other is Lagrangian formulation 

The first order in Eulerian formulation is well- 



y,|15|,|16|,|17|,y,|19|. 



21 



known linear theory and that in Lagrangian formulation corresponds to well-known Zel'dovich 
approximation p6|,^. 

Since current observations seem to point to i7 < 1 and/or A 7^ [^], it is important to 
develop perturbation theories in general Q and A. Perturbation theory in Eulerian space for 
Friedmann models (Throughout this letter, I mean models with arbitrary Q and A = by 
Friedmann models and models with arbitrary Q and A by Friedmann-Lemaitre models) are 
considered by several people |2^, ^ ^ . Lagrangian-space counterpart is considered 

by some authors [0,17,18,^, 2^, The second-order perturbation theories in Friedmann 



models can be expressed in closed forms. For Friedmann-Lemaitre models, however, the 
analytical expression for the time-dependence in second-order perturbation theories has not 
been known except for numerical solutions |30| , 15, 19, 22, 24 . 



In this Letter, I have obtained the explicit solution of the time-dependence in second-order 
perturbation theories both in Eulerian and Lagrangian space for the first time in models with 
arbitrary Q and A. 

Let us briefiy review second-order perturbation theories first. In Friedmann-Lemaitre mod- 
els, non-relativistic pressure-less self-gravitating fiuid are governed by the following continuity 
equation, Euler equation of motion and Poisson field equation 



5 + V ■[il + S)v] = 0, 

V + 2Hv + (v ■ V)v + V<l> = 0, 

2 



(1) 

(2) 
(3) 
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where v{x,t), $(£C,t) are position, peculiar velocity, peculiar potential in comoving coordi- 
nate, respectively, which correspond to ax, av{x, t), a'^^{x, t) in physical units, respectively. A 
dot denotes time derivative and V = d/dx denotes spatial derivative with respect to comoving 
coordinates. 

In Eulerian perturbation theory, these dynamical equations are solved perturbatively with 
respect to density contrast 6 and velocity v, keeping only terms of leading time- dependence 
The perturbative expansion is denoted as 5 = 5^^^ + 5^"^^ + ■ ■ ■, v = v^^^ + v^'^^ + ■ ■ ■ 
where (5'-"^ and v*^"^ are of order (^'^^^^ ■ The perturbative solution of the velocity field derived 
from this procedure is irrotational, V x v^^^ = 0, in the leading time-dependence of v^"\ 
Thus the perturbative velocity field is determined by its divergence 9^"''^ = V • v^^-yH as 
= HV A~^9^'^\ where A^^ denotes the inverse Laplacian. The result, up to second- 
order [|3],|3g,|l]^ 



IS 



6^'\x,t) = D{t)e{x), (4) 
9^'\x,t) = -f{t)D{t)e{x), (5) 
6^'\x,t) = D\t) + K{t)]e\x) + Ve{x) ■ V^{x) + ^[1 - ir(t)]v9,,(a^)^,,(a;)} , (6) 
0(2) = -f{t)D\t) [C(t)e\x) + Ve{x) ■ V^^^x) + [1 - C {t)]^ ,,{x)v ,,{x)] , (7) 

where e{x) corresponds to the initial fluctuation and ip{x) is its potential field, V'^ip = e and 
/ = D/{HD), K = F/D'^ C = F/{2DD) are defined using the growing solution D, F of the 
following ordinary differential equations: 

D + 2HD - -H^QD = 0, (8) 

F + 2HF - -H^QF = -H^QD\ (9) 
2 2 ^ ' 

The solution of the linear growth rate D governed by Eq. (8) is well-known. In expanding 
universes, it is 

D{t) oc an dxp/x + Ax^ + 1 - ^] - A]"^/^ (10) 
Jo 

Note that = 8ttGp/{3H^), A = A/{3H^) are time-dependent parameters. The function / is 
similarly expressed The solution of Eq. (9) for Friedmann models was explicitly given P^ . 
Although in fi-|-A = 1 models, this equation have been investigated numerically [^,^, I will 
give the closed representation of the solution of this equation for arbitrary fl and A universes 
below. 
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Lagrangian perturbation theory considers motion of mass elements labeled by unperturbed 
Lagrangian coordinates q. The comoving Eulerian position of mass element q at time t is 
denoted by x{q,t). The displacement field ^{q,t) defined by 

x{q,t) = q + ^{q,t), (11) 

is a dynamical variable in this formulation. Density contrast and a velocity field are derived 
from a displacement field as 

6[x{q,t),t] = J-'-l, (12) 
v[x{q,t),t]=x{q,t)=^{q,t), (13) 

where J = det[dxi/dqj] = det[6ij + d^i/dqj] is the Jacobian from Lagrangian space to Eu- 
lerian space. A dot denotes time derivative fixing Lagrangian coordinates q, so should not 
be confused with the Eulerian time derivative in Eq. (l)-(3). These Eulerian Eqs. (l)-(3) are 
transformed to the Lagrangian equations governing displacement field as 

JV^ ■ + 2H^] + ^H^VL{1 - J) = 0, (14) 
X [ip- + 2H^] = 0. (15) 

Again, a dot denotes Lagrangian time derivative and Vx is the spacial derivative with respect 
to Eulerian coordinates x, so {'Vx)i = J~^JijC^q)jj where Jij are cofactors of the Jacobian J. 

In usual treatment of Lagrangian perturbation theory, one assumes an additional condition, 
i.e., vorticity-free condition: 

Vx X = Vx X = 0. (16) 

This is not unreasonable requirement because vorticity is expected to be diluted by expansion 
before 'turn-around' in some sense P3[. The condition (16) is a sufficient condition for dynam- 
ical equation (15) [|T3[, so the solutions with vorticity-free condition are in a subclass of the 
general solutions (rotational perturbation is argued in the literature [|r^,|IU[). This subclass is 
achieved if one limit the initial condition to irrotational field because irrotational velocity field 
at one time remains irrotational at later times from Kelvin's circulation theorem. Therefore, 
Eqs. (14) and (16) are solved perturbatively for displacement field ^ = ^^^^ + ^^"^^ + • ■ •, 
keeping only terms of leading time-dependence. For irrotational initial displacement field, the 
result, up to second-order (see |jl9|,|3 for detail), is 

= -D{t)VA-h{x), (17) 
^(2) = J-K{t)VA-' E (vpffvI^W - v[;gvl.g) . (18) 
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The first-order solution is equivalent to the Zel'dovich approximation The factors 

D and K are the same quantities in Eulerian perturbation theory [Eqs. (4)-(7)]. Thus the 
importance of solving differential Eq. (9) in a closed form applies to second-order Lagrangian 
perturbation theory as well. 

The explicit form of the growing solution of differential equation (9) is now presented 
below. Bernardeau 0] realized that the solution of spherical collapse is useful in this problem. 
Although he claimed to obtain explicit results only for Friedmann models, I will show that 
this approach can lead to the explicit closed solution even in Friedmann-Lemaitre models. 

When the density contrast and the velocity field are spherically symmetric, Eqs. (l)-(3) 
or Eqs. (14) and (15) imply that the physical radius r of each mass-shell obeys the following 
equation: 



r = -—^ + - r, (19) 



GM A 
7^ ^ 3" 

where M is the total mass inside the mass-shell. This happens to be the same as Tolman- 
Bondi equation which describe the motion of mass-shell of spherically symmetric dust in 
general relativity, in spite of our Newtonian treatment. The density contrast at the center is 
evaluated by r ^ limit of the following quantity: 

2GM , , 

which satisfies 

A + - -H^QA = -H^A^ + (21) 

2 2 31 + A ^ ' 

I consider the solution A{t; 6i) with the initial condition, A{ti) = 6i and A{tj) = 0, which 
is expanded in Taylor series with respect to 6i: 

A{t; 6i) = D{t)6; + ^D2{t)6^ + ■■■. (22) 

The coefficient D obeys the differential equation (8) of linear growth rate and F = 31^2/2 — 21^^ 
obeys the differential equation (9). Thus, all we have to do is integrating Eq. (19) and then 
evaluating the coefficient of the Taylor series (22). The specific choice of the initial condition, 
A(ti) = 0, does not lose the generality because we are interested in a growing mode. 
Integrating Eq. (19), we have 

t-t, = j dr i^—— + -r^ + 2Ej , (23) 
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where a suffix i indicates initial values. E is an integration constant, 



2E 



2GM A 2 



l-VL-\ fi^i 



In the second line of the above equations, following equations are used: 



(24) 



f2GMY^\ 



-1/3 _ ^ 



Using a new integration variable x = {H'^^l/2GMy^^r, Eq. (23) reduces to 



Hit - U) = J 



/•(l+A)- 


1/3 






dx 








.X \ 



-2/3 



-1/2 



(25) 



Partially differentiating this equation with respect to 5i, the partial derivative dA/dSi is ob- 
tained as follows: 
dA 



— = (1 + Af/^^n + A(l + A)-i + A{1 + A)-V3 



X 



3M /-(i+A)-^/^^ 



5 - -— 



2 d6i J(l-H5i)-i/aai/a 



rfx ( — + Ax^ + A 

X 
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(26) 



where 



+ A(l + 5i)-i(ai/a)3 + A{1 + 5-:)-^l^a-Ja 

Partially differentiating Eq. (26) again, d'^A/dS^ is obtained similarly, though it is lengthy. 
Substituting Si — for these expressions, D and D2 are obtained explicitly. Note that Si — 
implies A = because of vanishing perturbation. 

In expanding universes, the term of leading time-dependence corresponds to the leading 
term of a/ai of the resulting expression. With this prescription, one can reproduce Eq. (10) 
for D as expected. Similarly, I find the exphcit expression for D2 and therefore for F. The 
final results are relatively simple: 



3/2 



3Ui 



1 - 



5/2 



2U: 



3/2 



c(n,A) = 



8/ 



-30 + 



2 + 4C/3/2-3(2 + n-2A)[/5/2]|, 



(27) 
(28) 
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where 

f/„(n,A)=/ dx\Vl/x + \x^ + (29) 
Jo 

These are the main resuhs of this Letter. In Einstein-de Sitter, Friedmann and flat models, 
f/„(l,0) = l/(l + a), ?7a(f^,0) =F(l,a,a + 2;l-^])/(l + a), and f/<,(f], 1 - fi) =F(l,a,(a + 
4)/3; 1 — i7)/(l + a), where F is the hypergeometric function. Moreover, in Friedmann models, 
f/3/2 and are are actually elementary functions which reproduce the previous results |16|. 
For practical purposes, Eqs. (27) and (28) are accurately approximated by 

K^^-n-^i^^-^^{i-\x\og,,n), (30) 

C^-^fi-"/^°°-4(l-^l°gio^)' (31) 

within maximum error 0.6% for both for —1 < log^g ^ < and < A < 1 . Quite obvi- 
ously from this fitting, dependence on parameters Vt and A is weak. In fact K and C for 
the above ranges of VL and A are different from the Einstein-de Sitter values within 8% and 
14%, respectively. Thus second-order effect is practically very insensitive either to or to 
A. Since this level of errors is negligible compared with the uncertainties with the current 
observational data, comparison of data and theories in a weakly nonlinear regime provides the 
test of gravitational instability almost independently of the poorly determined VL and A. Skew- 
ness [|,|16|,|6|,|3,||,|g,|3|,|32|, probability distribution function |g,|l|,|2[ and topology [||,|| 
in a weakly nonlinear regime are among such directly observable quantities. 

The analysis of weakly non-linear regime will be important as future redshift surveys will 
be available and statistical ambiguity on large scales will decrease. 
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